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Abstract 



We consider a general two-Higgs-doublet model with CP violation in the scalar 
sector. Three neutral Higgs fields of the model all mix and the resulting physical 
Higgs bosons have no definite CP properties. That leads, at the one-loop level 
of the perturbation expansion, to CP-violating form factors for jti, Zti and Wtb 
interaction vertices. We discuss asymmetries sensitive to CP violation induced by 

the form factors for the process e + e~ — > tt — > ■ ■ ■ and e + e~ — > tt — > b • • • at 
future linear e + e~ colliders. 
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1 Introduction 



Even though the top quark has been already discovered several years ago its in- 
teractions are still weakly constrained. It remains an open question if top-quark 
couplings obey the Standard Model (SM) scheme of the electroweak forces or there 
exists a contribution from physics beyond the SM. In particular, CP violation in 
the top-quark interactions has not been verified. The classical method for incor- 
porating CP violation into the SM is to make the Yukawa couplings of the Higgs 
boson to quarks explicitly complex, as built into the Kobayashi-Maskawa mixing 
matrix || proposed more than two decades ago. However, CP violation could 
equally well be partially or wholly due to other mechanisms. The possibility that 
CP violation derives largely from the Higgs sector itself is particularly appealing in 
the context of the observed baryon asymmetry, since its explanation requires more 
CP violation || then is provided by the SM. Even the simple two-Higgs-doublet 
model (2HDM) extension of the one-doublet SM Higgs sector provides a much 
richer framework for describing CP violation; in the 2HDM, spontaneous and/or 
explicit CP violation is possible in the scalar sector Q. The model, besides CP 
violation, offers many other appealing phenomena, for a review see Ref . |H . 

For our analysis, the most relevant part of the interaction Lagrangian takes the 
following form []: 

£ = ~ht{a + i-y 5 b)t + C^(m 2 z Z^ + 2m 2 w W^), (1) 

where h is the lowest mass scalar, g is the SU(2) coupling constant, v is the Higgs 
boson vacuum expectation value (with the normalization adopted here such that 
v = 2my/ / g = 246 GeV), a, b and C are real parameters which account for de- 
viations from the SM, a = 1, b = and C — 1 reproduce the SM Lagrangian. 

— CP — CP 

Since under CP, t(a + i^^,b)t — > t(a — i'-fsfyt and Z^Z^ — > Z^Z^ one can observe 

that terms in the cross section proportional to ab or bC would indicate CP vio- 
lation. The 2HDM is the minimal extension of the SM that provides non-zero ab 
and/or bC. 

In this paper we will focus on CP- violating contributions to the process e + e~ — > 

tt — > • • • and e + e" — > tt — ► b • • • induced within 2HDM. However the funda- 
mental goal is seeking for the ultimate theory of electroweak interactions. There 
are several reasons to utilize CP violation in the top physics while looking for 
physics beyond the SM: 

• The top quark decays immediately after being produced as its huge mass 
m t = 174.0 ± 3.2 ± 4.0 GeV JF| leads to a decay width T t much larger 
than Aqcd- Therefore the decay process is not contaminated by any frag- 

3 One could also consider more general, CP-violating ZZh coupling, see Ref. ||, however here 
the contribution from such a vertex would be negligible. 
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mentation effects [§] and decay products may provide useful information on 
top-quark properties. 



• Since the top quark is heavy, its Yukawa coupling is large and therefore its 
interactions could be sensitive to a Higgs sector of the electroweak theory. 

• At the same time, the TESLA collider design is supposed to offer an in- 
tegrated luminosity of the order of L = 500 fb -1 y -1 at */s = 500 GeV. 
Therefore expected number of ti events per year could reach 5 x 10 4 even for 
ti tagging efficiency e# = 15%. That should allow to study subtle properties 
of the top quark, which could e.g. lead to CP-sensitive asymmetries of the 
order of 5 x 10~ 3 . 

• Since the top quark is that heavy and the third family of quarks effectively 
decouples from the first two, any CP-violating observables within the SM are 
expected to be tiny, e.g.: i) non-zero electric dipole moment of fermions is 
generated at the three-loop approximation of the perturbation expansion , 
or ii) the decay rate asymmetry (being a one- loop effect) is strongly GIM 
suppressed reaching at most a value 10" 9 fllQf . So, one can expect that for 
CP-violating asymmetries any SM background could be safely neglected. 

Therefore it seems to be justified to look for CP-violating Higgs effects in the 
process of ti production and its subsequent decay at future linear e + e~ colliders. 
Even though 2HDM contributions to various CP-sensitive asymmetries has been 
already presented in the existing literature, see Refs. O, O, here we are providing 



a consistent treatment of CP violation both in the production, e e~ — ► tt, and in 
the top-quark decay, t — > bW. For an extensive review of CP violation in top-quark 
interactions see Ref. [p~3 



The paper is organized as follows. In Section |2|, we briefly outline the mechanism 
of CP violation in the 2HDM, introduce the mixing matrix for neutral scalars and 
derive necessary couplings. In section we present results for CP-violating form 
factors both for the ti production process and for t and £ decays. In Section |J we 
recall current experimental constraints relevant for the CP-violating observables 
considered in this paper. In Section we collect results for various energy and 
angular CP-violating asymmetries. Concluding remarks are given in Section [| 



2 The two-Higgs-doublet model with CP viola- 
tion 

The 2HDM of electroweak interactions contains two SU(2) Higgs doublets denoted 
by $i = 

{4>ti4 > V) an d ^2 — (02~; < / > 2)- It i s weu known [ffl that the model allows 
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both for spontaneous and explicit CP violation^. 

After SU(2)xU(l) gauge symmetry breaking, one combination of neutral Higgs 
fields, v^(c^30i + 3/3^2), becomes a would-be Goldstone boson which is absorbed 
while giving mass to the Z gauge boson. (Here, we use the notation sp = sin (3, 
= cos/5, where tan/5 = (0°)/ The same mixing angle, /3, also diagonalizes 
the mass matrix in the charged Higgs sector. If either explicit or spontaneous CP 
violation is present, the remaining three neutral degrees of freedom, 

<p 2 , ^ 3 ) = V2(W>i, ^2, spShfi - cpShfi) (2) 

are not mass eigenstates. The physical neutral Higgs bosons hi (i = 1,2,3) are 
obtained by an orthogonal transformation, h = Rip, where the rotation matrix is 
given in terms of three Euler angles (oti, a 2 , 013) by 

C\ — S\C 2 S\S 2 

R = ( S1C3 C1C2C3 - s 2 s 3 -C1S2C3 - c 2 s 3 I , (3) 
S1S3 C1C2S3 + s 2 c 3 -C1S2S3 + c 2 c 3 

where = sin a, and = cos ctj. 

As a result of the mixing between real and imaginary parts of neutral Higgs 
fields, the Yukawa interactions of the hi mass-eigenstates are not invariant under 
CP. They are given by: 

C = -^hJ(a{ + ib{ % )f (4) 

where the scalar (a{) and pseudoscalar (b{ ) couplings are functions of the mixing 
angles. For up-type quarks we have 

< = -R l2 , Vt = ^R t3 , (5) 

S/3 S/3 



and for down-type quarks: 



— Riii bf — —Ri3 , (6) 

C{3 Cp 



and similarly for charged leptons. For large tan/5, the couplings to down-type 
fermions are typically enhanced over the couplings to up-type fermions. 

In the following analysis we will also need the couplings of neutral Higgs and 
vector bosons, they are given by 



2 2 

HI HI 

g VVhi = = 2-^(spR l2 + c p R il ) 1 (7) 

for V = Z, W. Hereafter we shall denote the lightest Higgs boson by h and its 
R- matrix index by i. 

4 Here we are considering a model with discrete Z2 symmetry that prohibits flavor changing 
neutral currents. In order to allow for CP violation the symmetry has to be broken softly by the 
term /i^ 2 $}$2 in the potential. 
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Figure 1: Diagrams contributing to CP- violating form factors D^ z . 

3 Form Factors 
3.1 ti Production 

The effective tij and tiZ vertices will be parameterized by the following form 
factors^: 



Y(A V - B vl5 ) + {P \ Pl) \ c v - D vl5 ) 

ZUlt 



v(pt), 



(8) 



where g denotes the SU(2) gauge coupling constant, v — 7, Z, and 
4 . _ . v t a* 



A 1 = --sm9 w , B^ = 0, A z = — , 

3 2 cos W 

denote the SM contributions to the vertices for 



2 cos 6 



w 



C/D 



7,Z 



v t 



1 — - sin 2 9 W 



at 



The form factors A v , B v , C v describe CP-conserving while D v parameterizes CP- 
violating contributions. 

Further in this paper the following parameters will be adopted: m t = 175 GeV, 
m z = 91.187 GeV, T z = 2.49 GeV, sin 2 9 W = 0.23 and m b = 4.2 GeV. 

Since in this paper we are focusing on CP-violating asymmetries, the only 
relevant form factors are P 7 and D Z - Direct calculation of diagrams shown in Fig.|l] 



leads to the following result in terms of 3-point Passarino-Veltman [14] functions 



'Two other possible form factors do not contribute in the limit of zero electron mass. 
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defined in the appendix |A.l 



i Tn? 

9 

1 7TI 

D z = ^Az^bj[alm 2 t C 12 (pt,Pt,m 2 t ,ml,m 2 t ) 

-Cim 2 z C l2 {p u p h ml,m 2 t ,m 2 z )\. (9) 

Since the asymmetries we are going to discuss are generated by real parts of the 
above form factors, let us decompose (using eqs.(|5], 0)) 3?(-D 7 ) and $t(D z ) in the 
following wayQ 

3?£> 7 = R l2 R a f 2 t ] (10) 

3tD z = RizRizUm^ + /:23^ ^) + RilRi3f 13° ^ ) 

where superscripts indicate graphs that generate the contribution according to the 
notation of Fig.|I| Since in the case of the photon only Yukawa couplings a\ and 
b\ contribute, any signal of CP violation must be proportional to a\b\ ~ R i2 R i3 . 
However, for the Z boson vertex there is also other "source" of CP violation, 
namely b\Ci 3 RuR^. 

It is useful to discuss tan/3 dependence of the functions / first. From eq.(|]) 
and eqs.(^|, |7|) one can find out that all contributions to the form factors D 7 , 
Dz are enhanced for small tan/3: , ~ tan~ 2 /3, f 23 c ^ ~ cos/3 and 

fi3°~^ ~ tan _1 /3. Therefore, for tan/3 < 1, the contributions from diagrams c)-f) 
are expected to be suppressed relatively to those generated by diagrams a) or b) 
in Fig. [p. 

Hereafter we assume that there exists only one light Higgs boson h and possible 
effects of the heavier scalar degrees of freedom decouple. In Fig|J we illustrate 
dependence of the functions / on the lightest Higgs boson mass, m^. In order 
to amplify possible contributions we have chosen tan/3 = 0.5. As it is seen from 
the figure, the dominant CP- violating effects will be generated by and j 2 \ , 
which are generated by the Yukawa-coupling contribution, ajbj, from diagrams a) 
and b) in Fig.|I|. As one could have anticipated, we observe an enhancement of 
the contributions generated by the diagrams a) and b) for low Higgs boson mass. 
The growth of / 2 7 3 (a) and f 23 {b) is more pronounced for lower y/s (closer to the tt 
production threshold), it is a consequence of the Coulomb-like singularity generated 
by graph a) and b) in the limit rrih — > at the tt production threshold. Similar 
behavior have been also noted in the case of CP-conserving form factors |L5] . One 



also observes a typical threshold behavior at (y/s — m z ) for f 23 c ^ and f^ c ^ as a 
non-trivial absorptive part of diagrams c)-f) is necessary for nonzero contribution 

to gfcDsQ. 



The formulae for 5iZ? 7 and dlDz confirm results published in Ref . [[IT] . 



The same applies for 5RD 



7- 
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Figure 2: The functions / denned in eq. flTT| ) as a function of mh for tan f3 
y/s = 360 (left), y/s = 500 (middle) and yfs = 1000 GeV (right). 



0.5, 



If all three Higgs bosons of the model have the same mass, then the orthogonal- 
ity of the mixing matrix Rij guarantees [16| vanishing CP violation while summed 
over all the scalars. Therefore, the leading contribution to 9?-D 7 and $lDz originat- 
ing from an exchange of the lightest Higgs boson h could be partially cancelled by 
heavier scalars hh- However here we will assume that masses of the heavier scalars 
are above the production threshold for Zhh, therefore, as observed (for the lightest 
Higgs boson mass below 100 GeV) from Fig.|2|, the cancellation by heavier scalars 
for v/i = 360, , 500 and 1000 GeV could reach at most 20, 30 and 40%Q. 

Leading asymptotic formulae for small and large Higgs mass are presented in 
the appendix [A.2| . It is worth to notice that for non-zero f3 t = \Jl — Arn^/s, in the 
limit m\/ (PfS) — > both 3?D 7 and ^R.Dz are finite, whereas a typical decoupling 
limit, KZXp ^Dz ~ mf/(m\my / ) 1 is observed for large m^. 

Fig|| shows the functions / for tan/3 = 0.5 and two different Higgs boson 
masses = 10 GeV and rrih = 100 GeV, as a function of \/s. It is seen that it is 
not desirable to choose too high beam energy, as the size of the functions drops. 
Again, , dominate over , ^ by 1-2 orders of magnitude. 



8 Thus the cancellation expected in our case is not that strong as obtained in Ref. [jlj 
specific choice of model parameters. 



for a 
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Figure 3: The functions / defined in eq. 
m h = 10 GeV and m h = 100 GeV 



1 II j cLS db function of y/s for tan/3 = 0.5, 



3.2 Top Decay 

We will adopt the following parameterization of the Wtb vertex suitable for the t 
and t decays: 



iVi l Pl + f?P*) ~ '-^{^Pl + f 2 R P R ) 
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M w 
ia^k, 



(f 2 L P L + ffP R ) 



M. 



w 



111 



where Pl/r — (1 = F75)/2, Vt& is the (tb) element of the Kobayashi-Maskawa matrix 
and k is the momentum of W. In the SM f± = ff 1 = 1 and all the other form 
factors vanish. It turns out that in the limit of massless bottom quarks the only 
form factors that interfere with the SM are f% and f% for the top and anti-top 
decays, respectively. Currently, there is no relevant experimental bound on those 
form factors^. 

One can show that the CP- violating and CP-conserving parts of the form factors 
for t and t are not independent: 



fi' R = ±fi' U and 



L,R 



f L,R 
J2 



±f 2 R ' L , 



(12) 



where upper (lower) signs are those for CP-conserving (-violating) contributions |20], 
O . Therefore any CP- violating observable defined for the top-quark decay must 
e proportional to f l —fx or j 2 — / 2 ' . 



9 There exists direct experimental constraints from the Fermilab Tevatron on the form factors 
that are obtained through the determination of the VK-boson helicity. Pure V — A theory for 
massless bottom quarks predicts an absence of positive helicity W + bosons, therefore the upper 
limit on the helicity J-+ implies an upper limit on the V + A coupling f^, however, the resulting 
limit is rather weak pq| . There exist an indirect, but much stronger bound []l9| on the admixture 
of right-handed currents, ff 1 , coming from data for b — > 57, namely —0.05 < < 0.01. 
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Figure 4: Diagrams contributing to f 2 



Diagrams contributing to CP violation in the decay process are shown in Fig.^. 
Direct calculation leads to the following result for the CP-violating part of f 2 0: 



fR\ 

J 2 \cpv 



J 22 



^22 + C 



[2af6|C? 2 + (a\b\ - a$) CgJ - b\C t C 
;(^-l)+2(CS-CS)^ 
+ [2a\C*C{ 2 ct p - C f 22 Z(y!Cntp + C f 23 {Q(^C+)ct p + <Z{y\Cr)t p )\ 
~2a\C 2 C{ 2 ct^ - C 9 22 %y\C~)t p + Cf 3 (^Cf) $ct p + Q^C")^)] } 



(13) 



where 



C 



Cij(Pt,-Pw,m 2 h ,m^m 2 b ) 
CJp, 2 2 2 



- l3 \Pw,Pb,m h , m wi m t 
CJpw 



Pt,ml,m 2 H ,m 2 b ) 



Qij{pwiPhm h ,m H ,m t \ 



and 



t/b t/b , .jt/b 

Vi = a i + l W > 



where = spRi\ — cpRi2 and Q = i?^. The asymmetries that will be discussed 
here depends on real parts of f 2 and f 2 . It is easy to note from eq.([i~3|) that 
imaginary parts of C functions contribute to the real part of f 2 \cpv- It is seen 
that only diagrams b), d) and f) will contribute to %l{f 2 \cpv)- However, b — > 57 
strongly suggest [f?TJ that for 2HDM m#± > m 4 — m&, therefore eventually (adopting 
the relation (|i~2"|)) one gets the following result (from graphs b) and d) only) for 
CP-violating contribution to ^t{^j 2 \cpv)'- 

s g m b _ 



ta = tan p, eta = cot (3, 



(14) 



(15) 



CPV) 



167T 2 V 



22 



(16) 
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The general result for fi \cpv agrees with formulae for Qfi lcpv from Ref. |l 
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As will be discussed in Section ^ there is a strong experimental bound on \d\ for 
nih < 105 GeV. Taking into account the limit on \Ci\ and choosing tan/3 = 50 (in 
order to illustrate a possible enhancement) we plot dt(f^— f%) m Fig. [|as a function 
of nih- It is seen that ^(f^ — f^) is by 2—4 orders of magnitude below 9ft.D 7 or IfoDz 
even for large b-quark Yukawa coupling, compare Fig. ^ and [| The suppression 
is caused both by the experimental limit on |Cj| (for rrih < 105 GeV) and by an 
extra suppression factor of (mb/m t ) 2 (relative to ^D 1)Z )- 

There is a comment in order here; since in the 2HDM the real part of CP- 
violating form factors in the top decay is much smaller then in the production 
process, it is interesting to look closer at the suppression mechanism and find class 
of possible extensions of the SM that provide large 3?(/| ? — f%)- Since an absorptive 
part is needed, the only graphs that may contribute are those denoted by b), d) 
and f) in Fig.£| (allowing for the neutral scalar to be replaced by a neutral vector). 
One source of the suppression is the bottom-quark mass that originate from the 
propagator while the second one comes from the Yukawa vertex. The latter one 
could be easy amplified by large tan/?: for tan/5 ~ 38 the bottom-quark Yukawa 
coupling is as strong as the SU(2) gauge coupling. Therefore the suppression to 
overcome is m& from the bottom-quark propagator. A possible solution [22| seems 



to be a multi-doublet-Higgs model that could evade the stringent restriction from 
the b — > 57 decay and also overcome (through a contribution from the graph type 
f) in FigfJ) the limit on W + W~h that comes from the LEP limit on \Ci\. Let 



us notice that it is much easier to develop large ^(/g ~~ fk)i see e -S- Refs. |L2 



and [[T3[ . 

It is worth to mention that even though in the SM there exists one-loop contri- 
bution to 3?(/^ — /f), it turns out to be strongly GIM suppressed |TTJ. Therefore, 
although the 2HDM prediction for 3?(/^ — f%) is smaller than for CP violating 
from factors in the production mechanism, it is still by a factor ~ 10 5 larger then 
the SM result. 



4 Experimental Constraints 



Hereafter we will focus on Higgs boson masses in the region, rrih = 10 4- 100 GeV. 
As it has been shown in the literature [23] the existing LEP data are perfectly 
consistent with one light Higgs boson within the 2HDM. It turns out that even 
precision electroweak tests allow for light Higgs bosons E3[ . 

In order to amplify the form factors calculated in this paper we have adopted 
for an illustration tan/? = 0.5. However, there exist experimental constraints 
on tan /3 from K° — K° and — mixing [ , b — > s'j decay |^TJ and Z — > 
bb decay p6| . Since small tan/3 enhances H^tb coupling, in order to maintain 
tan/? = 0.5 we have to decouple charged Higgs effects and therefore we assume 
that m H ± > 500 4- 600 GeV. 

The constraints on the mixing angles a« that should be imposed in our numerical 
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Figure 5: 9ft (/^ — f%) as a function of rrih for tan/3 = 50. It has been assumed 
(according to constraints from b — *> 57) that the charged Higgs boson is heavy, 
m H ± > m t — rrib, therefore there is no absorptive part necessary to develop a non- 
zero contribution from diagram f) in Fig.|]. Additionally, we assumed Cj = Ri2 
(what is a very good approximation for tan/3 = 50) and selected such values for 
Ri2 and that are consistent both with Rij orthogonality and Cj LEP bound, 
and provide a maximal value of 3ft (Z^ — ZaO- 
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analysis are as follows: 



The ZZh couplings, Cf, are restricted by non-observation of Higgs-strahlung 
events at LEP1 and LEP2, see Ref. Wi 



• The contribution to the total Z- width from Z — »■ Z*hi —> ffhi is required to 
be below 7.1 MeV, see Ref. @. 

It turns out that the restriction on the ZZh coupling from its contribution to the 
total Z-width is always weaker then the one from Zh production if > 10 GeV. 

The LEP constraints on the ZZh coupling restrict the following entries of the 
mixing matrix Rif 

| sin /3i? i2 + cos^l <C- xp , (17) 

where C% xp stands for the upper limit for the relative strength of ZZh coupling 
determined experimentally in Ref. [27] up to the Higgs mass = 105 GeV. As we 



have concluded in the previous section, CP-violating phenomena we are considering 
are enhanced by small tan/?, in that case one can see from eq.(|T?J) that the LEP 
constraints mostly restrict Rn. Through the orthogonality the restriction on Rn 
is being transfered to constrain \R i2 Riz\ = |-R«2\A — RJi ~ ^ti\ which multiplies 
leading contributions to all CP-violating asymmetries considered here. The final 
result for upper limit on \Ri 2 Ra\ as a function of tan/5 is shown in Fig.||. In 
fact the bound on \R i2 Rjz\ depends on the Higgs mass, however, in order to be 
conservative, we have assumed C% xp = 0.12 that is the most restrictive experimental 
limit (obtained for m h ~ 18 GeVQ). 

As it is seen from Fig|] the constraints for \Ri2Ri3\ are weak for small tan/?. 
Therefore for tan/3 ~ 0.5 it should be legitimate to assume \Ri2R%s\ — 1/2 which 
is the maximal value consistent with orthogonality. 

Using the maximal value of Ri 2 Ra allowed by the orthogonality and the LEP 
constraints for small tan/5 = 0.5, we may discuss a possibility for an experimental 
determination of the calculated form factors at future e + e~ colliders. A detailed 
discussion of expected statistical uncertainties for a measurement of the form fac- 



tors has been performed in Ref. ||29|| . It has been shown that adjusting an opti- 
mal e + e~ beam polarizations, using the energy and angular double distribution 
of final leptons and fitting all 9 form factors leads to the following statistical er- 
rors for the determination of CP-violating form factors: A[3?(Z? 7 )] = 0.08 and 
A[3?(D^)] = 14.4 for e t i — 15%. It is seen that only 5J(.D 7 ), could be measured 
with a high precision. We have observed in Figs.|],|3| that 9?(-D 7 ) may reach at most 
a value of 0.10, therefore one shall conclude that several years of running with 

ii p or nih ~ GeV the limits presented in Fig. 16 of Ref. |27[] for the case when no 6-tagging 
and with 6-tagging almost coincide. Therefore our plot in FigB is not influenced by potential 
problems concerning the dependence of the Higgs-66 and Higgs-T + r~ branching ratios on the 
mixing angles. 
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10 1 10 tan/3 

Figure 6: Maximal value of \R i2 Ri^\ allowed by the LEP constraints on ZZhi 
coupling as a function of tan (3. 
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yearly integrated luminosity L = 500 fb _1 y _1 should allow for an observation of 
dt(D^) generated within 2HDM, provided the lightest Higgs boson mass is not too 
large. On the other hand, the expected precision for the determination of the 
decay form factors is much more promising: A[3?(/^ — f£)] = 0.014. However, as 
we have seen in Fig.|5], the maximal expected size of 3?(/^ — f£) is 5 x 10~ 5 (for 
rrih > lOGeV), therefore either an unrealistic growth of the luminosity, or other 
observables (besides the energy and angular double distribution of final leptons) 
are required in order to observe CP-violating from factors in the top-quark de- 
cay process. The results of Ref. J2I| assumed simultaneous^ determination of all 
9 form factors, therefore another chance to reduce of A[3?(/^ — f£)] is to have 
some extra independent constraints on the top-quark coupling coming from other 
colliders, like the Fermilab Tevatron or LHC. 



5 CP- Violating Asymmetries 



Looking for CP violation one can directly measure J29| all the form factors including 



those which are odd under CP. However another possible attitude is to construct 
certain asymmetries sensitive to CP violation. In this section we will discuss several 
asymmetries that could probe CP violation in the process e + e~ — > ti — > I ■ ■ ■ . We 
will systematically drop all contributions quadratic in non-standard form factors 
and calculate various asymmetries keeping only interference between the SM and 
D y , DzOTSUfX-ff). 



5.1 Lepton-Energy Asymmetry 

Let us introduce the rescaled lepton energy, x, by 



x = 



2E l fl-^ 1 ' 2 



m t \ 1 + p t 



where Ei is the energy of / in e + e cm. frame and j3 t = a/1 — Amf/s. Using lepton 
energy distribution da ± /dx calculated for the general form factors given in 



eqs. flgJlTD one can define the following energy asymmetry: 

., . . da~ I dx — da + 1 dx 
Acp{x) = - — -r~i , ' , 19 

CFy ' do~ I 'dx + da+ / 'dx v ' 

Direct calculation leads to the following result in terms of the CP-violating form 
factors: 

A i M _ 2^g(x) + m2 R -f2 L M(x) + r ] 5g(x)] 
ACP[ ) ~ 2[f(x)+r ]g (x)} • (20) 



2 Obviously, that leads to reduced precision for the determination of the form factors. 
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where 



for 



(3 - (3 2 )D V + 2(3 2 D A 



sin 9 W 

s 



5R 



(v 2 + a 2 )^ 



[s — m|) 2 64 sin 3 9 W cos 3 # 



ir 



Ulr 



V e V t 



I \ 16 sin 2 9 W cos 2 



6 sin cos 9 



-Di 



7 u 



£>v = (v e v t d - -) 2 + (a e v t df 
D A = (v e a t d) 2 + (a e a t d) 2 , 



with the SM neutral-current parameters of e and t: v e = —1 + 4 sin 2 9w, 
v t = 1 — (8/3) sin 2 ^^, and a 4 = 1, and a Z-propagator factor 

s — m| 16 sin 2 9 W cos 2 

The coefficient r\ is defined as 



V = 

for 



4^(1} 



(3 - (3 2 )D V + 2(3 2 D A 
2. 



D V a = v e a t d(v e v t d - ^) + (a e d) 2 a t v t . 



The definitions of the functions /, g, 5/ and <5g could be obtained from Ref. |30 . 

In order to estimate a relative strength of various sources of CP violation it is 
worth to decompose the asymmetry as follows: 

A 1 cp(x) = g\ ti {x) KD, + g l zti (x) ^D z + g l wtb (x) 3?(/ 2 R - f 2 L ). (21) 

As one can see from Fig.|7] the CP- violating effects that originate from the decay 
are substantially enhanced in the soft energy region. It is worth to notice that 
the minimal lepton energy for \fs = 500 GeV is E] nm ~ 7.5 GeV that is large 
enough to detect the lepton. Therefore the region of soft leptons should be carefully 
studied experimentally. The enhancement is a consequence of particular behavior 
of f{x) + rjg(x), g(x) and Sf(x) + r]8g(x) that causes the relative amplification 
of the decay effects, see Fig.l in Ref. []3l] . The same figure explains the observed 
smallness of the Z-boson contribution. For hard leptons both 3?(/^ — /jf) and 3?Z^ 7 
are enhanced and they are raising with the cm. energy. 

The energy-asymmetry could be decomposed into the leading contribution pro- 
portional to Ri2Ri3 and the remaining piece proportional to RuR^. The former 
one (that provides the leading contribution) is plotted in Fig.[8] for a fixed energy, 
x = 0.8, as a function of m^. One can observe that the largest asymmetry for the 
chosen energy corresponds to \fs = 360 GeV and m h = 10 GeV. 
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Figure 7: The coefficient functions g defined by eq.(^TJ) for the energy-asymmetry 
as function of x for yfs = 360 (left), 500 (middle) and 1000 GeV (right) for tan/3 = 
0.5. The solid curve represents the coefficient gL^(x), dashed g l zt ^(x) and dotted 

9wtb( x )- 
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Figure 8: The Higgs mass dependence of the coefficient of i?i2-Ri3 for the asymmetry 
given by eq.(0) for ^=360 (solid), 500 (dashed), 1000 GeV (dotted) for a fixed 
energy, x\ = 0.8, and taiaj3 = 0.5. 
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5.2 Integrated Lepton-Energy Asymmetry 

CP symmetry could be also tested using the following leptonic double energy dis- 
tribution 1021: 



' (l " r (22) 



a dx dx 



where x and x are for l + and I respectively, for 

d = l, c 2 = £, c 3 = ^(/ 2 fl -/^) 



and 



where 



fi{x,x) = f(x)f(x)+r]'g(x)g(x)+r}[f(x)g(x)+g(x)f(x)], 
f2(x, x) = f(x)g(x) - g(x)f(x), 

f 3 (x, x) = 6f(x)f(x) - f(x)5f(x) + rj'[6g(x)g(x) - g(x)5g(x)} 
+rj[Sf(x)g(x) - f(x)8g(x) + 5g(x)f(x) - g(x)6f(x)], 



, 1 {1 + /3 2 )D V + 2f3 2 D J 



(3 2 (3-(3 2 )D v + 2f3 2 D A ' 
The following asymmetry could be a measure of CP violation: 

dxdx-^— / / dxdx- 




tii ■/ ■/ .,- dxdx J J x >x dxdx /qq\ 

Acp = ' ' 3^ + / / ixii ■ (23) 




As before, it is useful to separate contributions from various form factors: 

A" P = <4 K-D 7 + g ll ztl KD Z + ^ K(/ 2 R - f 2 L ) . (24) 

In Table [l] we show the coefficients g for various cm. energies. Firstly, is clear 
that for any given a/s the coefficient g ll ztl is the smallest one. Secondly, it is seen 
that just above the threshold for tt production there is an enhancement of relative 
contributions from the decay, however that still not sufficient to overcome the 
suppression of 3?(/^ — f%) that we have observed in Fig.|5|. Therefore we can 
conclude that the leading contribution is provided by CP violation in the ^tt vertex. 



Fig.p illustrates the Higgs-mass dependence of the leading (proportional to 
R12R1?) contribution to the integrated lepton-energy asymmetry. It turns out that 
\fs = 500 GeV provides the largest asymmetry. 
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Figure 9: Higgs mass dependence of the coefficient of R^Riz for the asymmetry 
given by eq.(||) for ^=360 (solid), 500 (dashed), 1000 GeV (dotted) for tan/5 = 
0.5. 
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v^GeV] 


y-ytt 


a 11 - 


9wtb 


360 


0.0509 


0.00954 


0.410 


500 


0.386 


0.0684 


0.291 


1000 


0.602 


0.102 


0.235 



Table 1: The energy dependence of the coefficients g defined in eq.(|24|). 



Using results of Ref. 32 one can find out an expected statistical error for 
the determination of Aq P at any given e + e~ collider. Assuming y/s = 500 GeV, 
L = 500 fb y -1 an d lepton tagging efficiency, t\ = 60% we get AA" P = 0.014. As 
it is seen from Fig.^| an observation of the asymmetry would require several years 
of running at the assumed luminosity. 

5.3 Angular Asymmetry 

Another CP- violating asymmetry could be constructed using the angular distribu- 
tions of the bottom quarks or leptons originating from the top-quark decay: 



da 



d cos 9) 



= 3n ^™ Bf (n f + n{ cos e f + n f 2 cos 2 o f ) , (25) 



where / = b,l, Bf is an appropriate top-quark branching ratio, 9f is the angle 
between the e~ beam direction and the direction of / momentum in the e + e~ cm. 
frame and fl{ are coefficients calculable in terms of the form factors, see Ref. [B3 
The following asymmetry provides a signal of CP violation: 



da + %) da-(TT-6 J 



1) 



^U0f) = tTl /r 9 ' „ (26) 
cpyu da+(9 f ) | da-{ir-9 s ) > 

d cos 6 f d cos Of 

where da + /~ is referring to / and / distributions, respectively. Since Of — >• ix — Qj 
under CP, the asymmetry defined above is a true measure of CP violation. 



Adopting general formulas for the asymmetry from Ref. [S3| and inserting form 
factors calculated here we plot the asymmetry in Figs.pl], |IT| as a function of cos 
for bottom quarks and leptons, respectively. As before, the asymmetry can be 
decomposed into ^ti, Zti and Wtb vertex contributions: 

Mjp%) = 9U°f) ^ + 9 f ztM) ^D z + g f wtb (0 f ) U(f« - f 2 L ). (27) 

It is seen that forward-backward directions are favored, however an experimental 
cut | cos#/| < 0.9 should be imposed in the realistic experimental environment. 
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1 COSl^ 




Figure 10: The coefficients functions g denned by eq.(f27|) for the angular asym- 
metry, for final bottom quarks (/ = b) as a function of cos 6 for ^/s = 360 (left), 
500 (middle) and 1000 GeV (right) for tan/3 = 0.5. The solid curve represents the 
coefficient g^ tl (6 f ), dashed g Zt ti^) an< ^ dotted g^^iOf) 
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Figure 11: The same as Fig.llO] calculated for final leptons (f — I). 



In order to illustrate the Higgs mass dependence we plot in Fig.|12| the angular 
asymmetry both for f = b and / = I for chosen polar angle cosQ = —0.25. 
As expected the maximal effect could be reached for minimal Higgs mass rrih = 
10 GeV, \fs = 500 GeV is the most suitable energy. 



5.4 Integrated Angular Asymmetry 

The angular distribution given in eq.(^) could be adopted to define an integrated 
[29H of the angular asymmetry A{j P {O f ): 



version 



/ a cos v f — — a cos 9, 

J~c m dcosO f J 



d cos Of 



A f (P- P +) = m } ±2 1 (o c 

/ dcosOf— ^+ / dcosOf— 

J dcos0 f J dcosOf 

where P e - and P e + are the polarizations of e and e beams, da + ^~^ is referring 
to / and / distributions respectively, and c m expresses the experimental polar- 
angle cut. In order to discuss possible advantages of polarized initial beams we are 
considering here dependence of the asymmetry on the polarization. Hereafter we 
will discuss the same polarization for e and e: P = P e - = P e +. 
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Figure 12: The Higgs mass dependence of the coefficient of Ri2R^ for the angular 
asymmetry defined by eq.fl2~6|) for v /i=360 (solid), 500 (dashed), 1000 (dotted) 
GeV for a fixed polar angle cos# = —0.25 and tan/3 = 0.5. 
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V^[GeV] 


P 


quark b 


lepton 








9ztt( p ) 


9wtb(P) 




9ztt( p ) 


360 


0.0 


0.00844 


0.00106 


0.142 


-0.0162 


-0.00203 




0.8 


0.00983 


-0.00555 


-0.259 


-0.0493 


0.0278 




-0.8 


0.00758 


0.00510 


0.388 


-0.0106 


-0.00713 


500 


0.0 


0.113 


0.0136 


0.121 


-0.224 


-0.0270 




0.8 


0.131 


-0.0718 


-0.247 


-0.627 


0.343 




-0.8 


0.101 


0.0661 


0.347 


-0.149 


-0.0968 


1000 


0.0 


0.332 


0.0389 


0.0678 


-0.722 


-0.0845 




0.8 


0.422 


-0.225 


-0.167 


-1.55 


0.824 




-0.8 


0.284 


0.181 


0.194 


-0.507 


-0.322 



Table 2: The energy and polarization dependence of the coefficients g^ a (P), g za {P) 
and gly tb {P) defined in eq.(|^) for leptons (/ = I) and bottom quarks (/ = b). 



Again we decompose the asymmetry as follows: 



^cp(p) = gU p ) ^ + 9L(P) ^Dz + gU(P) - (29) 

In Table ||] we show the coefficient functions g calculated for various energy and 
polarization choices assuming the polar angle cut | cos 9 f \ < 0.9, i.e. c m = 0.9 in 
eq.(p8|), both for leptons and bottom quarks^ . It could be seen that a positive 
polarization leads to higher coefficients <r tf and g Zi ^. Since 9ft(-D 7 ) > 3t,(Dz) ^> 
9R(/:f — fz) that implies that maximal asymmetry could be reached for P = +0.8 
and the dominant contribution is originating from 3ft(Z? 7 ). Since the number of 
events does not drop drastically when going from unpolarized beams to P = +0.8, 
it turns out that the positive polarization is the most suitable for testing the 
integrated angular asymmetry. It is clear from the table that the asymmetry for 
final leptons should be larger by a factor 3 + 4 then the one for bottom quarks and 
their signs should be reversed. 



Using the general formula for the asymmetry from Ref. |29] and adopting results 
for the CP- violating form factors we plot ^4^, p (P e -, P e +) in Fig.|T3| as a function of 
the Higgs mass both for bottom quarks and leptons. It is clear that the largest 
asymmetry could be expected for P e - = P e + = +0.8 for final leptons at y/s = 
500 GeV. With the maximal mixing, R i2 Riz = 1/2 the 1% asymmetry could be 
expected for the Higgs boson with mass = 10 + 20 GeV. Since the statistical 
error expected p9| for the asymmetry is of the order of 5 x 10~ 3 , we can conclude 



that the asymmetry A f CP (P e -, P e +) is the most promising one, leading to 2a effect 

13 Note that in Tabic || there is no column corresponding to the coefficient of 3?(/^ — /^). 
That happens since the angular distribution for leptons is not influenced by corrections to the 
top-quark decay vertex, see Refs. |33L M and 129], 
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Figure 13: The Higgs mass dependence of the coefficient of i?j2-R«3 for the angular 
asymmetry defined by eq.(p8|) for bottom quarks (upper) and leptons (lower) at 
^=360 (solid), 500 (dashed), 1000 GeV (dotted) with unpolarized beams (left), 
P = +0.8 (middle) and P = -0.8 (right) for tan/3 = 0.5. 



for light Higgs mass and tan/3 
have polarized e + e~ beams. 



0.5. As it is seen form Fig.|13| it is relevant to 



6 Summary and Conclusions 

We have considered a general two-Higgs-doublet model with CP violation in the 
scalar sector. Mixing of the three neutral Higgs fields of the model leads to CP- 
violating Yukawa couplings of the physical Higgs bosons. CP-asymmetric form 
factors generated at the one-loop level of perturbation theory has been calculated 
within the model. Although in general the existing experimental data from LEP1 
and LEP2 constraint the mixing angles of the three neutral Higgs fields, their 
combination relevant for CP violation is not bounded for small tan f3 which is the 
region of our interest. We have shown that the decay form factors are typically 
smaller then the production ones by 2-3 orders of magnitude. The dominant con- 
tribution to CP violation in the production is coming from ytt coupling. Several 
energy and angular CP-violating asymmetries for the process e + e~ — > tt — > ■ ■ ■ 

has been considered using the form factors calculated 



and e + e 



tt 



► b 
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Figure 14: Convention for momenta and mass labelling used. 

within the two-Higgs-doublet model. It turned out that the best test of CP invari- 
ance would be provided by the integrated angular asymmetry A f CP (P e ~, P e +) for 
positive polarizations of e + e~ beams. For one year of running at TESLA collider 
with the integrated luminosity L = 500 fb~ y -1 one could expect 2a effect for the 
asymmetry for light Higgs boson and tan/? = 0.5. 
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A Appendix 

A.l 1-Loop Integrals 

Here we will define 3-point one-loop integrals used in calculations of form fac- 



tors. The convention for momenta and mass labelling is presented in Fig. 14. The 



Passarino-Veltman functions [|14[] C are defined as follows: 



// 



4— n 



d n k 



1, k^, k^k v 



(27r) n [k 2 - m\\ [{k + Pi) 2 - m| [{k + pi+ p 2 ) 2 - mi] 

1 Q IJZ „2 „2. _2 _2 _2a 



16tt 2 



I A I I I I I 

{p 1 ,p 2 ,p ;m 1 ,m 2 ,m 3 



(30) 
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The vector and tensor integrals and C can be expanded into scalar coefficients 
and Lorentz covariants: 

C = KC n +P^C 12 (31) 

cr = p^c 21 +p^c 22 + (p^+pMc 23 + ^c 24 . (32) 

Imaginary part of C 12 for certain sets of arguments is needed for the calculation 
of K£> 7 and MD Z : 



7T 



^C 12 (p t ,pt;m t ,m h ,m t ) = — 

spt 



-tt log 1 + 



" ft 2 



7T 



%C 12 (p t ,Pi;m h ,m t ,mz) = x 
[2t 2 (l - ft 2 - , 2 ) - ft 2 /3 2 ] I log * - ^ 



ft 2 + A/3z 



0(^- 2771*) 

(33) 

- P z > ©(v^ -m h - m z ), 



for ft = nth/ y/s, p t — \/l — 4m 2 / s and /3z = A(l, z 2 , ft 2 ), where A(a, 6, c) is the 
standard kinematic function. 

In order to calculate Re(f 2 R — f 2 ) one needs an imaginary part of the following 
C 22 function. Neglecting m b and defining m w = and m h = for uo 
1 — rn^jm 2 one gets: 



uimt 



3C 22 (p w ,p b , m 2 w , m 2 h , m\ 



71LU 



2mf 



{ [2m 2 (6rh 2 w + 1) - (2m 2 w + 1)] + 



-2 [m 2 h (6m 2 w + 1) + 2m 2 w ] m\ log (l + mf) } . 



(34) 



A. 2 Asymptotic Formulae 
A. 2.1 Production of ti 

For ft, = ^ in the limit ^ < 1 one gets: 



t 2 m 



27T/3 t V 2 



, ft 2 , 

1 V og 



0(1 - 4t 2 



(35) 



1 m 



2ir$ v 2 



'AzbAatfp 2 



1 + ^ 



A 2 



9(1 -4t 2 )+ 



- C^ 2 



A(l - ft 2 - ^ 2 ) + (2t 2 (l - ft 2 - z 2 ) - h 2 p 2 t ) log 



1-Pt 



0(1 -ft -2),} (36) 



where z = m z j \fs and t = m t j \J~s. 
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For -§r ^> 1 one gets: 

p t s 



s ^ = -i4^ , "4( 3 -4) e ( i -« 2 » < 37 > 



A. 2. 2 Decay of the Top Quark 

In the limit m h = < 1: 



ujmt 

W?\cpv) = ^ (^) 2 m W { -(1 + +[6 + m-J + 2 log(m^)] } , 



while for ^> 1: 



/ ys 9 ~^ \ / ^ ^ 2 1 , A o 

-("V + 3) + (g"V + 



(38) 



(39) 
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